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Abstract 



We develop a relativistic covariant and unitary description of the pion-nucleon 
interaction in a hadron-exchange model. The model is based on the solution 
of a dimensionally reduced (quasipotential) Bethe-Salpeter equation for the 
partial-wave off-shell vrA^ scattering amplitudes with the potential consisting 
of the field-theoretical s- and it-channel nucleon, Roper, Delta, -D13, 5ii ex- 
changes, and the f-channel p and a meson exchanges. The contributions of the 
spin-3/2 Delta and Dis resonances are treated within the Rarita-Schwinger 
formalism and different forms of the vrA^A vertex are investigated. The free 
parameters of the model are fitted to the vrA^ phase-shift data of the KH80 
and SM95 partial-wave analyses in the region up to 600 MeV pion kinetic 
energy. The resulting on-shell solution provides a good description of the vrA^ 
scattering lengths, as well as the energy behavior of the S, P, and D partial 
waves. The sensitivity of the phase shifts on various model-dependent effects 
is examined. 
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I. INTRODUCTION 



The TiN interaction received much attention in the past both theoretically and exper- 
imentally in view of its fundamental nature (early literature can be found in Refs. |[l|J^). 
Current theoretical interest is triggered by the experimental programs being carried out at 
NIKHEF, MAMI, TJNAF and other intermediate-energy facilities with the purpose to un- 
derstand the structure of hadrons and their interaction in the confinement region of QCD. 
To extract most physics from the new high-precision measurements a reliable and accurate 
knowledge of vrA^ and A^A^ interaction is required. Highly successful attempts have been 
made in describing these interactions in terms of hadronic degrees of freedom over a wide 
energy region. In particular, the relativistic one-ho son- exchange models were successfully 
applied to the description of A^A^ interaction, and especially during the past decade this 
theoretical framework has been extended to the vrA^ system [^-[TOl . Such an extension gives 
one a capability to study a broad class of reactions, including pion scattering and produc- 
tion on light nuclei in a self-consistent framework. In this paper we report on a relativistic 
covariant t^N model using the quasipotential approach and based on an effective interaction 
characterized by a hadron-exchange potential (some of our results have already been briefly 
reported 

Although the underlying dynamics of the vrA^ interaction is nowadays believed to be 
governed by QCD, it is practically impossible to resolve it fully in terms of quarks and 
gluons because of the confinement problem. Much of the present understanding of the vrA^ 
physics at low and intermediate energies remains to be based on dispersion relations and 
effective chiral Lagrangians [0,^ in term of the hadron degrees of freedom. 

The chiral pion-nucleon Lagrangians are usually extended in two ways: first, by includ- 
ing higher-mass states, such as p-meson, A-isobar, etc.; secondly, by including the higher- 
derivative terms. Both ways are necessary to extend such a phenomenological description 
to higher energies: The contributions due to higher-mass states have a clear physical sig- 
nificance, while the higher-derivative terms are needed to examine the effect of unknown 
short-range physics. The higher-derivative terms play, for instance, a crucial role in the 
renormalization program of chiral perturbation theory (ChPT). In the hadron-exchange 
models, a similar role is played by the "strong form factors" which are included in the ef- 
fective Lagrangian to model the short distance behavior of the potential. Both ChPT and 
hadron-exchange models thus begin from a similar "extended" chiral Lagrangian but the 
approach to calculating the vrA^ scattering amplitude is somewhat different. In ChPT one 
usually performs a perturbative field-theoretic calculation maintaining crossing-symmetry 
and exact agreement with the soft-pion theorems. (For the development of ChPT in ap- 
plication to the ttN scattering see Ref. [T^.) In the hadron-exchange approach one uses 
the effective Lagrangian to construct the potential which then is resumed via a scattering 
equation. In this way crossing symmetry is given up in favor of exact unitarity in a given 
channel space, and possibility of studying nonperturbative phenomena such as dynamical 
resonances. 

In defining a hadron-exchange model one usually specifies three ingredients: effective 
Lagrangian, potential, and the scattering equation. These ingredients are interrelated in 
quantum field theory, where one must solve the Bethe-Salpeter (BS) equation and has a 
well-determined procedure for computing its kernel from a given Lagrangian. The BS kernel 
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consists of all the irreducible graphs and hence can be computed only perturbatively. In 
this work we shall take the potential to be given by the tree-level graphs, although we 
do not claim that the used Lagrangian justify any perturbative expansion. On the other 
hand, the resulting approximation transparently relates to the usual quantum-mechanical 
picture where the scattering problem is given by a Lippmann-Schwinger type of equation 
for one-particle exchange potentials. Therefore, one might prefer to view this approach as 
relativistic quantum-mechanical one rather than some "nonsystematic" truncation of QFT. 

The four-dimensional BS equation for the vrA^ system (with a one-particle-exchange po- 
tential) has been solved by Nieland and Tjon |T^, and recently by Lahiff and Afnan |^ in a 
more realistic setup. Models P,^,|7|JT0| exploit instead various quasipotential (QP) equations, 
which can be obtained by a three-dimensional reduction of the BS equation. The use of QP 
equations provides a technical simplification of the problem, without destroying the Lorentz 
invariance of the theory. It should be remarked, however, that some of the QP equations can 
violate charge conjugation symmetry, and because this symmetry is crucial for renormalizing 
the positive- and negative-energy baryon poles in an equivalent way, the equations which 
violate it are less preferable. We will employ the equal-time (ET) equation which preserves 
the full Lorentz covariance, including charge conjugation. This equation will be specified in 
the next section. 

Apart from the technical simplifications, the QP approach can sometimes be motivated 
by physical arguments. For instance, while the four- dimensional BS equation for t-channel 
(meson-exchange) potential, i.e., the ladder BS equation, has a wrong one-body limit, a 



number of QP equations with the proper limit can be devised [|16H T9[|. In the vrA^ situation 



the potential may, in addition to the t-channel meson exchanges, contain the -u-channel 



baryon exchanges, which spoil the standard one-body limit arguments [2C]. The ordinary 
ET equation is shown to provide an optimal choice in the case when both t- and w-channel 
exchanges are present in the force corresponding to the vrA^ situation |[20|1 . 

As for the effective Lagrangian, we use the pseudovector nNN coupling, and include a, p 
mesons, and all the relevant (for the considered energy region) nucleon resonances as explicit 
degrees of freedom. The precise form of the Lagrangian and the potential is discussed in 
Sec. IV. An interesting aspect, which comes in with the nonperturbative modelling, is that 
for a sufficiently attractive potential the nucleon resonances can be generated dynamically. 



as quasi-bound states of the vrA^ system. In models [p!5|,pT| and quantum hadrodynamics 



2^ the A(1232) is described in this way. Lahiff and Afnan ^ include the A explicitly, but 



suggest that the Roper resonance can be of dynamical origin. Gross and Surya [^] include 
the A and the Roper poles, but treat the resonance dynamically. In this paper we con- 
sider P33A(1232), Pn A^*(1450), P>i3 iV*(1520), and 5iiA^*(1535) resonances. Within our 
model, these resonances are all of nondynamical origin, i.e., are included explicitly via an 
effective Lagrangian description. Of course, the dynamical effects will anyhow contribute 
to the generation of the resonances seen in the phase-shifts. Thus, an admixture of both 
"elementary" and "composite" component constitutes the full result. Since the "elemen- 
tary" fields corresponding to the resonances are included with real masses, the dynamical 
contributions are fully responsible for generating the width. Our model maintains the elastic 
vrA^ unitarity and therefore only the one-pion decay width of the resonances is generated. 

As we have to deal with the spin-3/2 fields of resonances, such as that of the A-isobar, we 
shall address here some of the problems of consistent formulation for relativistic higher-spin 
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fields. Consistent formulations for the free spin-3/2 field have of course been known for a 
long time. The Rarita-Schwinger formalism based on the vector-spinor field representa- 
tion became the most popular one. The form of the free spin-3/2 action is uniquely (up to 
trivial field redefinitions) constrained by requirements of Poincare invariance and consistent 
degrees-of-freedom counting. The latter requirement essentially means that the action must 
have enough symmetries to kill off the unphysical lower-spin components, and maintain only 
the physical 2s + 1 degrees-of-freedom of the theory. An arbitrarily constructed interaction 
may violate this consistency and activate the unphysical degrees of freedom. This neces- 
sarily leads to a number of pathologies, such as acausal propagations pi|,pB|, inadmissible 



quantization p6| , p7| , violation of Poincare invariance. 

The conventional vrA^A coupling |]28| , p9| ,3l is an example of such pathological interactions. 
The contribution of the unphysical spin-1/2 components appear in the scattering amplitudes 
via the dependence on the so-called "off-shell parameter" and "spin-1/2 backgrounds." On 
the other hand, a class of consistent vrA^A and 7A^A couplings has recently been established 
30| , |3ll| . Those are essentially all possible couplings that maintain the gauge symmetry of the 



free massless Rarita-Schwinger action. In the present model we shall use the leading "gauge- 
invariant" vriVA coupling, which has the same non-relativistic limit as the conventional one. 
In Sec. |V| we study the differences between the conventional and the gauge-invariant ttA^A 
coupling at the tree-level A-exchange contributions. The largest differences are seen first 
of all in the spin-1/2 partial waves, where the conventional coupling gives the background 
contributions verse no contribution from the consistent couplings. 

The parameters of the effective Lagrangian, including the form factor masses, form the 
set of model parameters. Some of them, such as the irNN coupling constant, the nucleon, 
and the meson masses, are very well-determined elsewhere and therefore are kept fixed 
during the fits. The others are fitted to give the best agreement with the vrA^ partial-wave 
analyses ||3^-p4| . The complete model provides an accurate description of the S- and P-wave 
scattering lengths (x^/datapoint ~ 1.4), as well as the energy behavior of the 5*-, P-, and 
some of the D-wave phase-shifts up to 600 MeV pion lab kinetic energy. 

In the next section we describe the covariant quasipotential equation for the off-shell vrA^ 
amplitudes, and its partial-wave decomposition. This will specify the equation solved in the 
model. In Sec. |T| we briefly discuss the effective Lagrangian used to read off the tree-level 
potential — the driving force of the equation. The renormalization procedure to treat the s- 
channel poles is described in Sec. |^ In Sec. 0we analyze the effects of the different exchange 
contributions in the low-energy vrA^ data, and then present the results of the complete model 
emphasizing the effect of the rescattering contributions. Some discussion and concluding 
remarks are given in Sec. Finally, various appendices contain some technical details of 
the analysis. In Appendix A we summarize the conventions. Appendix B shows some details 
of the partial-wave and isospin decomposition of the off-shell vrA^ amplitudes. Appendix C 
provides explicit expressions for various hadron-exchange contributions to the off-shell ttN 
potential. Amplitudes for higher-spin baryon exchange are discussed in Appendix D. 
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II. QUASIPOTENTIAL APPROACH 



The fully off-shell relativistic vriV scattering amplitude in the space of the nucleon helicity 
spinors is described by sixteen scalar amplitudes: one for each combination of the helicity and 
p-spin of the initial and final nucleon. Parity conservation reduces the number of independent 
scalar amplitudes to eight. As a suitable covariant representation which expresses the off- 
shell amplitude in terms of eight invariants we choose the following: 



Tj:,,V,fc';p,fc) = <'(p)(i,j^') 



All Ai2 
A21 A22 



+ 



Bii B12 
B21 B22 



where Aij and Bij are the eight scalar functions of invariants formed by the momenta, i.e., 
p^, p -p', etc. The four-momenta of the initial (final) nucleon and pion are given by p and k 
{p' and k'), respectively, while P is the conserved total four- momentum of the system: 



P 



p 



k 



p 



k'. 



Due to the momentum conservation only three of the external momenta are independent, 
below we usually work with p, p' and P. Furthermore, are the nucleon helicity spinors, 
where A = ±i and p = ±1 (A' and p') are the initial (final) helicity and p-spin of the nucleon, 
respectively. 



For the on-shell situation (p 



p 



m%, k^ 



k" 



ml 



the amplitude reduces to the 



standard form [p5|] : 



Tva(p',p;P) = 4'^(p') A{s,t) + \^-{k' + k)B{s,t) u'^'ip), 



(2) 



[s = P^,t= {p-p') 
and m 



A2 



U 



{p — k') are the Mandelstam invariants]. We find from Eqs. (|ip 



B = Bii+mNiBi2 + B2i) + m%B22, (3) 
A = tunB + Aii+mNiAi2 + A21) +m%A22. (4) 

Our starting point for the vrA^ amplitude is the Bethe-Salpeter (BS) equation, schemat- 
ically shown in Fig. |l], 

•d^g 



T{p', p) = V{p\ P)+^J^ V{p', q) G{q) T{q, p) , 



(5) 



where V is the potential, G is the TtN propagator; the dependence on the total momentum 
P is omitted. If q is the relative four-momentum of the intermediate vrA^ state, the ttN 
propagator takes the following form: 



G{q) 



{(3P - g)2 -ml + ie {aP + qf 



{aP + g) ■ 7 ruN 

771% + is ' 



where 



a = a{s) = p ■ P/s = (s + 7n\ — m^)/2s, 
(3 = (3{s) = k-P/s={s-ml + ml)/2s. 



(6) 



(7) 
(8) 
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In approximating the BS equation one often simplifies the singularity structure of the 
kernel VG, such that the temporal integration can easily be done. This procedure is called 
three-dimensional (3D) reduction while the resulting equation is a quasipotential (QP) equa- 
tion. For instance, in the reduction to the spectator equation all the poles of V and 
the negative-energy pole of G in the go plane are neglected. 

As we have recently emphasized [llO| , p]T|| , the danger of doing a 3D reduction via ap- 
proximating the pole structure is that the charge conjugation symmetry can be destroyed. 
In particular, in our naive interpretation of the spectator equation this symmetry is vio- 
lated, essentially because of an asymmetric treatment of the positive- and negative-energy 
states. Gross has recently presented an interpretation of the spectator prescription which is 
consistent with the charge conjugation symmetry ||36|| . 

The equal-time (ET) reductions (see, e.g. such as Salpeter's instantaneous approx- 



imation |^8|, preserve charge conjugation symmetry. In these reductions one effectively 
removes the go poles from the potential while treating exactly the poles of the two-particle 
propagator G. To remove the potential poles one fixes the relative-energy variable go in 
some way. Most frequently the constraint go = 0, or its covariantized form: P ■ g = 0, is 
used. 

We will be using the ET type of approach. To implement the constraint P ■ g = 0, we 
may impose the condition that the interaction is insensitive to the off-shellness along the 
direction defined by an unit four- vector n^. For the two-body case this means that V and T 
entering the scattering equation depend on the projections of the relative four-vectors onto 
a 3D hyperplane orthogonal to n^. Defining the projection operator, 

Oi,y = g^u - n^n^, (9) 

we write the corresponding equation as follows: 

T(r, !) = V{1 + ^/0 Vil', q) G{q) T{q, I), (10) 

where /, g are the relative momenta of the initial, final, and intermediate irN state, re- 
spectively; = Ofj_J'^, and similarly for g. 

Equation ([T0| ) is manifestly covariant. On the other hand, it can easily be reduced to 



the 3D form. For instance, we can choose the frame where n = (1, 0, 0, 0), and therefore V 
and T are independent of the 0-th component of relative momenta (since any scalar product 
will depend only on the spatial components, e.g., g ■ 7 = — g ■ 7). The integration over go in 
Eq. ([T0| ) can now be readily done leading to the 3D equation. To prevent the dependence of 
the S'-matrix on n, one may choose n along some physical four-momentum, for instance along 
the total momentum: ra^ = P^/ V^P^. Then the reduction is possible in the center-of-mass 
system (CMS) where P = (Po, 0,0,0). 

In the ET approach, the two-particle propagator is sometimes modified to include ap- 
proximately the crossed graphs [p!8| , |39| , thus providing the correct one-body limit of the 
equation in the case of t-channel type of potential. We however do not apply such modifica- 
tions here, because they actually worsen the predictions for the vrA^ case where the w-channel 
exchanges are present |^ . 



Because of rotational invariance and parity conservation it is convenient to partial-wave 
decompose Eqs. (^ and (p!0| ) (see Appendix B). Let 
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n'-'xip',?) = Y^jiJ + l)Dixi^P'p)Tyi'ipip',Po,p;Po), (11) 

where Qp'p is the sohd angle between jf and p. Furthemore, in the CMS, using Eq. (|l]) and 
the Dirac equation 

- ptun) u1{p) = {po - pEp)-fo uKp), (12) 
the off-shell amplitudes can be written as 

n:',ip',p;Po)=4iP) [l+Ti''{p',p;P)+^^T^^{p',p;Po)] u{{p), (13) 

where 7± = |(/ ± 70), and T±'' are eight scalar amplitudes with definite parity ±. 

In doing so we in particular find for the case of the ET equations, that the parity- 
conserving amplitudes T/ satisfy 

00 

T^P'P = V.'P'P + -[dqq^T G^S V.'P'P" T/""", (14) 

where 

GS(|g1;Po)=2^ J ^ G^P\q,q,; P,) 
J zn 

—00 

= -p{ujg{-pPo + Eg + ujg-ie)y^ . (15) 

In this work we will be focusing on solving this equation for a one-particle-exchange vriV 
potential, see Fig. ^, where the potential is regulated by form factors. In the next section 
we describe the vrA^ interaction used in this study. 



III. EFFECTIVE LAGRANGIAN AND THE ttN POTENTIAL 

In the following we specify the interaction Lagrangian of vr, p, a, N and isobar fields, 
used to construct the hadron-exchange force depicted in Fig. ^ and written out in Appendix 
C. The field representation and the corresponding free Lagrangian is chosen according to 
the spin and isospin of the particle. Thus, the pion is described by scalar isovector multiplet 
7r°- = (tt"'', 7r°, 7r~), the sigma-meson is a scalar isoscalar field a, the A(1232) is represented 
by a vector-spinor isoquadruplet A^j = (A^"*", A^, A°, A~), and so on. 



A. Nucleon and meson exchanges 



The nNN interaction Lagrangian is taken in accordance with the chirally-symmetric 
In Weinberg's nonlinear realization the nN scattering amplitude to the 



a models 12,13 



leading order is given by the nucleon Born term with the pseudovector coupling plus the 



Weinberg- Tomozawa contact term |41,42|. The pseudovector coupling reads 
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4S = ^ iV7,75r.iV 9,7r^ (16) 

where /ttatat is the pseudovector vrA^A^ couphng constant (the pseudoscalar couphng constant: 
QttNN = '^f-KNN ij^N /f^-n) will also be used below). The Weinberg- Tomozawa contact term 
can be represented as a p-meson exchange with the following interaction: 

C,MN = 9pNN Nj {l,p'a + P^^j N, (17) 

^pnw = QpTT-K £abc Pa'^b 9^11^ (18) 



provided the p coupling, = gpNNgfm-K-, is fixed by the KSRF relation [^]: Qp = nip/ (v^/tt), 
where f^^ ~ 93 MeV. There is also the second form of the KSRF relation ||28|: Qp = 
V^mpfirNN/ (fn-KgA), dA — 1.26, obtained from the first one by using the Goldberger- 
Triemann relation for It should be remarked that the Weinberg- Tomozawa contact term 
is equivalent to the p-exchange only at threshold and provided gp is fixed by KSRF relation 
while Kp = 0. The energy dependence is different, but not significantly in the considered 
energy region. 

Since we use the pseudovector uNN coupling, the a-meson is in principle not needed 
from the standpoint of chiral symmetry. Nevertheless, a a exchange can be used to model the 
isoscalar contribution of the correlated two-pion exchange. In order to keep the agreement 
with the soft-pion theorem, a derivative coupling to the pion is used, i.e., 

jCaNN = QaNN (tNN, (19) 
Cann = CT d^TTi d'^Tli , (20) 

where the sign of the (TTrvr-coupling is chosen in accordance with the correlated two-pion 
exchange analysis 0, and is different from the one used in 0. This interaction leads to the 
following on-shell potential: 

y _ gaNNgamr t — 2m^ 

To control the effect of the cr-exchange on the scattering length we introduce a free 
parameter in the following way: 



Svrm^ m'i — t 



For Co- = 1 the a contribution to the S'-wave scattering lengths vanishes. Note that this 
modification amounts to adding the following term to the Lagrangian: 

^CTTrn = -TT CaTTl^a TT^ . (23) 
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B. A-isobar exchange and higher resonances 



The coupling of the spin-3/2 A field to the pion and the nucleon is conventionally de- 
scribed by the following Lagrangian (see, e.g., Ref. 0): 



C.NA = ^ \g^.u + i)7M7J Ta N S.vr'^ + H.c, (24) 



where z is the off-shell parameter, Ta is the isospin-i —>■ | transition operator. 

As remarked in the Introduction, this coupling involves the spin-1/2 sector of the Rarita- 
Schwinger field, and give rise to an unphysical spin-1/2 background. The latter effect can 
in principle be removed by inserting the spin-3/2 projection operator "by hand" in either 
the vertex or the propagator. For instance. Gross and Surya |Q have chosen this option. 
However, because of the nonlocal nature of the projection operators their use is problematic: 
unphysical singularities occur at s = and u = for the s- and w-channel contribution, 
respectively. In Ref. @ this problem is actually not met because the s = point is well 
below the threshold, while the w-graph vanishes in the approximation of that model. If the 
M-channel A exchange is present, one usually prefers to keep the background and fit the 
off-shell parameter . 



We shall also study the following ttA^A coupling pU|,pT| 



4na = ^^^'''^"^(5/.A.)757«T,iV9^7r« + H.c., (25) 

referred to as gauge-invariant ( GI) nNA coupling. Being invariant under the Rarita- 
Schwinger gauge transformation: 

A^(x) A^(a;) + d^e{x), 

where e{x) is a spinor field, this coupling does not involve the spin-1/2 components of the A 
field. As a consequence, the spin-1/2 backgrounds are totally absent from the corresponding 
A-exchange amplitudes. 

We include also the s- and w-channel graphs of Pn (Roper), ^n, and Dy^ resonances. 
At low energies the contribution of these resonances is marginal,^ but they are important 
for the proper description at higher energies. The first two particles are treated same as the 
nucleon but with different masses, coupling parameters, and, in the case of Su, parity. The 
Di3 is treated in the same way as the A (the same propagator and interaction vertex), but 
with different isospin, parity and mass. Exchanges of even higher spin resonances can in 
principle be easily included in our model via the amplitude obtained in Appendix D. 



^This is generally not true for a spin-3/2 resonance if the conventional coupling is used, since the 
spin-1/2 background can be large even far away from the mass position. 
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C. Cutoff form factors 



The high-energy behavior of the hadron-exchange field-theoretic potential is usually reg- 
ulated using the off-shell form factors introduced in the vertices. We have introduced them 
for each of the particle in the vertex. For the pion we use the monopole form factor: 

= (26) 

For the a- and p-meson we use the one-boson-exchange form factor: 

Uit) = -7^- (27) 



For the baryons we use the form factor of Ref. 



Mr) = "„f J (28) 



with n = 2. 

In addition, for each pion we introduce the following cutoff: 



W(^'^) = 7r^- (29) 



This function is motivated by considering the effect of the higher-mass states on the high- 
energy behavior of the ttN propagator G. If we were to include not only the pion and 
the nucleon but all the states lying on the infinitely rising Regge trajectories, at very high 



energies their effect factorizes in the form of function like Eq. (^) The energy depen- 
dence and the absence of singularities distinguish this cutoff from the usual ones, such as 
the monopole form. We take the same value A^, for the and f Regge cutoff masses. These 
two form factors do not affect the on-shell potential, because no pion exchanges appear in 
the Born graphs. 

It is important to realize that the final results depend on the off-shell form factors, 
even after the renormalization is applied. The physical meaning of such form factors is 
usually given in analogy with that of the electromagnetic form factors. They thus reflect the 
extension of the hadrons, and in principle should be calculated from the underlying theory. 

Our flt to the vrA^-scattering phases will determine the values of the cutoff masses. They 
are given in Table [V| together with the rest of the model parameters. Using these values, in 
Fig. ^we have plotted the form factors which affect the loop contributions. Their dependence 
on the loop momentum is shown, while the 0-th component is flxed by the equal-time 
constraint and the energy is flxed at threshold. 

The actual cutoff of the model is given by the solid line in the flgure. As one can see, 
it is rather soft: it starts off as a monopole with the mass about 0.8 GeV, and is even 
softer above = 0.5 GeV^. At higher energy it becomes softer as well, because f Regge is 
energy-dependent. However, the latter effect is small as can be seen from Fig. ^, where the 
energy dependence of f Regge is shown over the region of our vrA^ flt. The energy dependence 
of /Ar(s) and /a(s) is shown there as well. 
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IV. RENORMALIZATION 

Since there are s-channel singularities in the considered potential, we have to carry out 
a renormalization procedure. We adopt the scheme in which the Lagrangian is expressed in 
terms of the physical parameters and no "bare" parameters appear. Then, in principle, the 
counter-terms should be subtracted and fixed by the renormalization conditions. To perform 
such a renormalization procedure it is convenient to work with the one-particle-irreducible 
Green functions. One can separate the potential into two terms V = Vg + Vu, where 

Vs{p',p) = T.W) SB{P)rB{p) (30) 

B 

represents the s-channel baryon exchanges (pole terms), while Vu contains the rest of the 
graphs (nonpole terms). Since VgG is a separable kernel, we can explicitly resum these 
contributions, and find that the resulting amplitude can equivalently be written as 

T{p\p) = Y: rUp') S*BB'{P) nip) + Tu{p',p), (31) 

BB' 

where 

nip) = nip) + ^ /0 n{q) G{q) Uq,p), 

i^BB') ^ = iS^)^^SBB' - ^BB', 

Sbb' = I j'^,rB{k)G{k) r^,{q), 
and Tu satisfies the following integral equation, 

Tu{p\p) = Vu{p\p)+i I^Vu{p',q)Giq)Uq,p). (35) 

The full amplitude is thus written in terms of the irreducible Green functions: S*,r* 
and Ty. The diagrammatic form of this representation is given in Fig. ^. 

A. Baryon mixing 

Note that the dressed baryon propagator, Eq. (^3]), is in general nondiagonal. In other 
words, the baryons can mix. Of course this mixing happens only among the baryons with 
the same "good" quantum numbers, such as spin and isospin. Parity is also conserved, 
nevertheless the mixing of baryons with the same spin, isospin and opposite parity may 
occur due to the negative-energy state propagation. 

To perform the renormalization we need first to diagonalize the propagator. Since it is 
a complex symmetric matrix, we diagonalize it using a complex orthogonal transformation 
O [pO^ = O^O =1). The full solution can obviously be written in the diagonal form as 
follows: 

T = Y^ (ro)^ {o'^s*o)^ {o'^r)^ + r,, (36) 

B 
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(32) 
(33) 
(34) 



In our model we include only two baryons with the same spin and isospin (nucleon and 
A^*). For this case the propagator is diagonalized by a 2x2 complex orthogonal matrix which 
can be parametrized as usual by one complex variable, 

0=( ^^"^V (37) 

Y — sm X cos X J 

where in this way we introduce the A^A^* mixing angle x- 

Furthermore, since we use the same Feynman rules for the nucleon and A^*, their dressed 
vertices are equal up to the coupling constants. Therefore for the A^A^* self-energy matrix 
one can write 



/ ^7VAf(-P) ^NN*iP) ) _ ( dwNN gnNNdnNN 

while for the vertex 



E{P), (38) 



"^^1 = 1 3ttNN 



r. (39) 

The propagator is then diagonalized by the orthogonal transformation (R^) with 



X(P) = \ arctg ( 2 f ^ - ^ - "^"'"7,^. ) " | . (40) 

I \9ttNN* QttNN g-nNNQ-nNN" ^\J^ j / \ 

The corresponding eigenvalues are clearly 

Sn{P) = [f- rriN - Sn{P) - ^e]-' , (41) 
Sr,*iP) = [f - m^* - Sn^P) - , (42) 



where 



Sn = {gnNN COS X + girNN* siu xf ^ + {^11^* - mn) SVD? X, 
^N* = {gnNN* COS X " gnNN siu xf ^ ~ i^N* - niN) SIV? X- (43) 



The vertices are rotated according to 



Pn = {gnNN COS X + g^NN* siu x) P* , 

Pn* = {g^NN* COS X - g^NN sin x) P*- (44) 



B. Self-energy 

Let us consider the mass renormalization using the counter-term method. The counter 
terms can be read off directly from the free Lagrangian. For the spin-1/2 case, for instance, 
they are given by Z2{mQ — m) -|- (1 — Z2){f— m), where mo is the bare mass, and Z2 is the 
field renormalization constant. The renormalized spin-1/2 baryon propagator is defined as 
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ie\ 



(45) 



where r^^"(f) = E{f) - ^2(7710 - m) - (1 - m). 

In the CMS, P = {Pq, 0), the self-energy can be written as 

r-"(Po) = i^r(Po)7+ + i^:^"(Po)7-, 

and a similar decomposition holds for the propagator, 

S(Po) = SW(Po)7+ + S(-)(^o)7- 

where 7± = (J ± 70), and 



(46) 



(47) 



SW(Po) 
S(-)(Po) 



Po - m - r;'^"(Po) - is 
Po + m + ri^"(Po)+ie 



(48) 



Obviously, 7+ and 7_ act as the projection operators onto the positive and negative energy- 
states, hence 8^+^ corresponds to the positive and S^~^ to the negative energy-state propa- 
gations. 

The renormalization condition at the pole position is given by 



(Po-m)S(+)(Po; 
(Po + m)S(-)(Po: 



Po=m 
Po=—m 



Expanding Z'+(Po) near Pq = m, and Z'_(Po) near Pq = 
condition requires 



(49) 



-m, we find that the renormalization 



dS+{Po 



dPo 



Po=m 



dPo 



Po=-m 



As emphasized earlier [|l^ , the above described procedure breaks down if the self-energy is 
computed using a quasipotential formulation which violates charge conjugation symmetry, 
since in that case Z'+(Po) 7^ Z'„(— Pq). The self-energy of the spin-3/2 baryons can be 
renormalized similarly, since the spin-3/2 baryon contribution to the spin-3/2 partial- waves 
can always be factorized into vertices and a spin-1/2 propagator, see Eq. ( P7|) . 



C. The renormalized vertex and the amplitudes 

In the adopted renormalization scheme we require (i) the (real part of) renormalized 
baryon self-energy Z'^^"(Po) and its first derivative vanish at the pole position Pq = ±m; 
(ii) the (real part of) renormalized vertex ttN ^baryon vertex is equal to the undressed 
vertex at the renormalization scale n defined as the point where all three particles are on 
the mass-shell, fi: k"^ = m^, = m^, P^ = m^. 
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For the vertex we use the multiphcative renormahzation since it maintains unitarity in 
a simple way. The renormahzed vertex is thus defined as 

Tip;P) = Z^rip;P), (50) 

where Zi is the couphng constant renormahzation factor which is readily determined from 
condition (ii): 

z, = r(/i)/r*(^). (51) 

In the case of the A^A^*-mixing we renormalize the (scalar) function S in Eq. (^) at the 
point associated with the nucleon. This procedure clearly yields the proper physical nucleon 
mass pole in the corresponding baryon propagator. Adopting this subtraction procedure the 
A^* mass and coupling constant at the nucleon mass position can then be extracted. In the 
various tables the values of these parameters found in the fits are quoted. 

After the partial-wave decomposition, the renormahzed solution of the ET equation for 
a given isospin I and total spin J and parity r reads as follows (for brevity the external 
momenta are omitted): 

T/"'^ = J2 S^'^"^^ 6j^j 6j^j + Ti^;p, (52) 

OO 

p" 

n = ( r/' + i Jdgg^ E rf{q) G^^'^iq) Tl^'^iq)] , (54) 
\ ^ p' / 

Sr = z,- [dqq'J2rnq)G^^kq)niq), (55) 

where Jb, Ib and tjb are the baryon spin, isospin and parity, respectively. The renormahzed 
propagator is 

SW(Po) = [±Po -rn~ r-"(Po) + te]~' , (56) 

where the renormahzed self-energy is given in terms of S of Eq. (|55D [with Eq. (^) in the 
case of A^A^*-mixing] as follows: 



(57) 

Po=±m 



V. RESULTS 

Having described the equation for the off-shell vrA^ amplitudes, its renormahzation, and 
the driving force, we now turn to discussing the outcome of such modeling. For this let us first 
give an explicit relation between the on-shell amplitudes + + and the phase parameters. 
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We introduce the standard on-shell amplitudes fi±^ where / = J — |r is the angular 
momentum and r is the parity of the state. In the normalization according to 

where 5i± is the phase shift and rji± is the inelasticity, we can identify 

fi± = q '^""l"""' T/ + -^(E, g, g; Po), (59) 

where 

q = ^J[s- {tun - m^)2] [s - {tun + m^)2]/4s, (60) 



and E = \Jm\ + is the energy of the nucleon in the CMS. 

At very low energies the partial-wave amplitudes are dominated by the threshold behav- 
ior: g^^"*"^, and their real and imaginary parts are related by elastic unitarity. Therefore, it 
is sometimes more useful to study |^| 



/ \ 21+1 

Mum = [y] (^^^ 

instead of fi± itself or the phase shifts. Note that, at g^ = 0, 7Vl/j is equal to the corre- 
sponding scattering length defined as0 

au^hm[q-''-'fi±m. (62) 
The nN effective-range parameters hu can also be determined in terms of ^A■. 

»u-^,Mu(e%__^. (63) 

Instead of using this formula, we will be presenting the plot of 7W as a function of g^. 
The slope of these "effective-range plots" at small g indicates the values for the effective- 
range parameters. In the partial waves that support a resonance (e.g., Pn, P33) it is more 
appropriate to study another effective-range expansion: 

g2'+i cot 5i± = — + \rij g2 + ■ • • , (64) 

nontheless, we only will address threshold parameters a and h for all partial waves. 



^We shall commonly refer to this quantity as to scattering length, even though for P and higher 
waves it is properly called scattering volume because of the dimension. 
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A. The /C-matrix approximation 



In this subsection we focus on the K-matrix approximation to the full scattering problem. 
Results of the full model are discussed in the next subsection. In the K-matrix approximation 
vrA^ amplitude is given by the lowest order i^'-matrix expression: 

where Ki± = qa{s) V^, / = J =F 5, and l^"' is the partial-wave potential [obtained from the 
potential by means of Eqs. (|1^) and ( pi2|) ]. Equation (|65D clearly satisfies elastic unitarity, 
but the principal value of the loop integrals is neglected. 

This approximation is considered to be good, at least at low energy, and has been fre- 



quently used, see e.g. ||4^ , ^ for most recent applications to the vrA^ scattering. At low 
energies, indeed, the soft-pion theorems dictate that the Born graphs dominate, implying 
for the potential modelling that the rescattering effects should be relatively small. When the 
latter is true, considering the K-matrix approximation may allow us to make a preliminary 
adjustment of some model parameters without going to the full calculations. We in partic- 
ular would like to examine the effect of using different ttNA couplings in the A exchange 
contribution. 

As is well known, the S-wave scattering lengths are well reproduced by the Born-level 
nucleon and p- meson exchanges alone [^,0. Taking 91^^/4:71 = 13.6 and (7^/4?? = 3.0, we 



plot in Fig. y the contribution of these three graphs (s- plus u-channel nucleon exchange 
plus t-channel p exchange) to Aiu for all the 5* and P partial waves. The figure shows that 
the ^-wave scattering lengths are indeed reproduced. However, the energy dependence of 
the 5* waves is not well described. In fact, the slopes of the effective-range plots have a 
wrong sign. Furthermore, the P waves are not reproduced at all. Similar picture occurs if 
the contact term is used instead of the p exchange. 

Looking at the ratio of the S-wave lengths it is certainly plausible that they should be 
dominated by some isovector contribution,^ such as the p-meson exchange. Therefore, it 
would be interesting to find a simple mechanism which accounts for both the P waves and 
the energy behavior of the 5* waves, and, at the same time, does not affect the S-wave 
scattering length. Since P33 has the largest discrepancy, we study first the effect of the 
A-isobar exchange. 

In Fig. 1^ we show the calculations performed with the two different choices of the off-shell 
parameter: 

z = —\ (Peccei choice [^), 
z = \ (NEK choice [||), 

and /^tva/^tt = 0.36. The nucleon and p couphng constants are kept the same as in the 
previous calculation, with = 3.7. 



^ Experimentally the ratio 0511/0531 ~ —1.75, which is close to the ratio of the isospin factors for 
an isovector meson contribution. 
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One can see that, as far as P33 is concerned, the A contribution is very plausible for both 
choices of z. The large difference between the choices clearly shows up in the spin-1/2 partial 
waves, where the A exchange produces a significant background contribution controlled by 
the off-shell parameter. 

In particular, Peccei's choice affects substantially the S waves, and hence spoils the 
scenario where those are dominated by the p-exchange. The P waves look much better, 
which could be due to the remarkable fact that, for this particular choice, the s-channel 
A-exchange graph gives only a tiny contribution to the spin-1/2 P waves. The NEK choice, 
in contrast, does not affect the S waves at threshold and gives a large effect in the P waves. 

Clearly, Peccei's choice could be favorable phenomenologically as long as the missing 
strength in the S waves is somehow explained; for instance, by an isoscalar meson exchange. 
We believe such a scenario is realized in most of the models which use coupling p^ ) and 
describe the scattering lengths correctly. However, apparently it is not possible to describe 
simultaneously the S- and P-wave scattering lengths in the tree-level model with only the 
p, and A exchanges. 

The Gl ttNA coupling Eq. (|25|), in combination with the usual Rarita-Schwinger prop- 
agator, leads to the A-exchange amplitude with only spin-3/2 contributions (omitting the 
isosp in-dependent factor) : 

"—ill (4) 7^^"'^* 

where pu = P ~ k — k' , and 

is the spin-3/2 projection operator. 

A calculation using this amplitude is shown in Fig. |^, in comparison with analogous 
calculations using the conventional coupling with Peccei and NEK choices of the off-shell 
parameter. First of all we remark that the A contribution to the spin-1/2 partial- waves 
comes from the -u-graph only, and not from the spin-1/2 components. From the figure we can 
see that the GI and NEK coupling produce similar contributions to the S'-waves, but largely 
different results in the spin-1/2 P- waves. The resonant P33 wave comes out very much alike 
for both GI and conventional couplings, the main difference being {s/m\) factor in front of 
the GI amplitude. Hence, at threshold the GI result is factor of {rriN +m^)^/m^ ~ 0.76 
smaller than the conventional result (modulo small contributions from the w-channel graph, 
which for instance are responsible for the difference between the Peccei and NEK choice in 
the P33 wave). Despite that, after a readjustment of parameters the GI invariant coupling 
usually gives a better description of the scattering lengths, see calculations presented by 
Tables | and |I| (note that here the NEK choice has been used with Hp = 0, as is suggested 
in the original paper and indeed gives a better fit than with Kp = 3.7). The problem 
with the wrong energy behavior of the S waves, however, applies to all these calculations. 
To correct for this a scalar cr-meson exchange is needed in the ttN force. 
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Inclusion of the a-exchange allows us to fit the scattering lengths to practically arbitrary 
accuracy, independently of whether we use the conventional (model A) or the GI coupling 
(model B), see Tables fT| and |V]. Since we have fixed = 1, and thus do not allow the 
(J to affect the S'-waves, the best fit of the off-shell parameter give the NEK value which 
also has vanishing S-wave contribution. The S-wave lengths are therefore explained in both 
models A and B in exactly the same way: by p and nucleon exchanges alone. The small 
difference between the NEK and GI coupling in P33 (see Fig. 0) is apparently compensated 
by the difference in Jt^na- The large differences in the other three partial- waves has mostly 
been removed by taking a different value for Kp. This example shows that if the potential 
is general enough, the spin-1/2 backgrounds of the A can possibly be reshuffled into other 
contributions. 



B. Full calculations 



We solve the ET equation, Eq. ([T^), by Fade approximants following the procedure 



described in Refs. Writing the equation as 

T = V + VGT, (68) 

we begin by performing several iterations of the potential, and hence find the first few terms 
in the expansion of the amplitude: 

00 00 
T = Y, ^^"^ = V+J2 ViGV)''. (69) 

n=0 n=l 

The solution is then sought in the form of the Fade approximant. For the equal-time equation 
with the model potential the solution accurately converges by performing just six iterations. 



We have fitted the on-shell solution to the KH80 [g and SM95 fS^ nN scattering 
partial-wave analyses, in the region from the threshold up to 600 MeV pion kinetic energy 
in the lab. The resulting fit is shown by the solid lines in Fig. ^. The determined coupling 
constants and masses are given in Table |V]. 

The dashed line in the Sn phase-shift of Fig. ^ indicates the calculation without the 
s-channel Su resonance graph. This graph contributes also to the Pn wave but calculation 
with or without it produce practically identical results. The Su resonance pole is thus 
relevant only for the Su wave above 400 MeV. 

Up to 350 MeV the agreement of the model with the partial-wave data is very good as 
can be seen in Fig. ^. At energies exceeding 350 MeV inelastic channels become important. 
Since we have not considered any inelastic mechanisms, some discrepancies seen in Fig. § at 
higher energies are not surprising. 

In fitting we have paid careful attention to the low-energy behavior, particularly to 
the correct shape of the "effective-range plot", A^/±(g^), defined in Eq. (|6lD. The model 
description of these plots is shown in Fig. |TU[ The scattering lengths can be read off these 
plots at g = 0. The x-square value for the scattering lengths with respect to the SM95 
analysis is 1.4. 

To give a feeling about the size of the rescatterings in the model, the dashed lines in 
the figures indicate the results of the calculation where the principal part of the rescattering 
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integrals is neglected, i.e., the i^-matrix approximation. Unlike in the ii'-matrix calculations 
of the previous subsection, here the form factors are included, and the same set of parameters 
is used as in the full calculation. 

The large difference between the full and the i^-matrix calculation in the waves where the 
baryon pole contributions are present indicates, therefore, that the "dressing" and renormal- 
ization of the pole contributions has a significant effect. The effect of the rescatterings in the 
nonresonant waves, such as 5*31, P31, etc., which have contribution only from the T„ term, 
is smaller. We also observed that the "attractive waves" (which have positive scattering 
length) receive comparably large and positive rescattering contributions at threshold. 

For the resonant waves the T„ contribution may lead to significant shifts of the resonance 
position. For instance, the A-pole in the P33 is located at ^/s = = 1.252 GeV, while 
inclusion of leads to the physical P33 wave which resonates at y/s = 1.232 GeV. In the 
i^-matrix approximation the resonance always occurs at the pole position. To describe P33 
in this approximation one is to use ttia =1.232 GeV. 

By readjusting the parameters we are able to reproduce the phase-shifts in the i^'-matrix 



approximation up to 350 MeV, see Fig. 11. The values of the parameters are given in 



Table 0. Note, however, that using the full model we could achieve the fit of a better 
quality and up to higher energies. 

The (s-channel) baryon pole terms are modified by the vertex and self-energy corrections. 
As has been seen from comparing the full and i^-matrix calculation, such a "dressing" may 



have appreciable effects, especially in the resonance Pn and P33 waves. In Fig. |1^ we plot 
the real and imaginary part of the nucleon and the A isobar self-energies. From the figure 
we can see that the energy dependence is indeed significant. The same is observed for the 
NN* mixing angle plotted in Fig. 0. The renormalization of the pole terms produces the 
values of the renormalization constants given in Table |V14 

The vertex corrections are studied using the dynamical form factor defined as follows: 

FP{q^s)=TP{q',s)/rP{q',s), (70) 

where F'' is the undressed off-shell vertex, and F'' is the renormalized off-shell vertex, see 
Eq. (0). Note that the coupling constants and the cutoff form factors are cancelling out in 
the expression ([70D, since they are the same for both of the vertices. The dynamical form 
factors are thus equal to unity at the renormalization point. 

The model prediction for ttNN, ttNA and nNDi^ form factors is given in Fig. O as a 



function of the off-shell 3- momentum for ^/s = mN + m.,^, and in Fig. as the function of 
energy for the on-shell situation, \q\ = q. According to these figures the rescatterings have 
much larger effect on the energy dependence then on the off-shell behavior of the vrA^ state. 
Recently, the uNN form factor has been studied by Saito and Afnan PS|, and by Schiiltz, 



Haidenbauer and Holinde in a similar modelling. In the latter work the ttA^A form factor 
has been studied as well. To compare with the results presented in we need to multiply 
our dynamical form factor by the off-shell form factor given in Fig. |^ (solid line). We then 
can see that the resulting irNN form factor of our model agrees in the main features with 
that of model EStl. It is therefore less soft than the form factor found in Ref. EHI- This 



allows the rescattering contributions to play a bigger role. 
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VI. DISCUSSION AND CONCLUSION 



Throughout the calculation we have been fixing the ttNN coupling constant to the value 
advised by the Nijmegen group [Q: fl^Nl^'^ — 0.0757. Our fits were not very sensitive to 
the increase of the coupling towards more traditional value of 0.078. 

The value of Qp comes out close to the one inferred by the p-meson decay width: gp/4:7r ~ 
2.8. It is also consistent with the KSFR relation (Sec. IV A), which gives 2.78 in its first 
form, and 2.94 in the second form. The small value of Hp has been mainly dictated by the 
simultaneous fit to S^i and P31 waves at the higher energy scale. At low energies the change 
in Kp affects mostly the Pn channel, as can be seen from Fig. ITDI where the dash-dotted line 
represents the calculation with the vector meson dominance value: Kp = 3.7. 

Comparing Tables |V| and|VJwe see that, depending on whether the rescattering contribu- 
tions are included or not, very different values of the A-isobar masses and coupling strengths 
are needed to obtain correct phase shifts. This indicates how the dynamical component due 
to ttN loops may play a significant role in the generation of the observed A(1232) resonance. 
That is in addition to the "elementary" component due to the formation of the three-quark 
state. 

In comparing with other relativistic models we can comment that our major difference 
with the model of Pearce and Jennings |^ resides in the A-isobar and cr-meson contributions. 
For the A, they use the conventional coupling, and for the crvrvr coupling they use Eq. (pOD 
with the opposite sign and without the additional Ca- term Eq. (P^D- The a contribution 
is thus attractive in their case and can account for the discrepancy in the S-waves which 
appears due to the A background. Although it probably gives rise to problems at higher 
energies, which is indicated by the very soft cutoff form factor for the a, with cutoff mass 

= 0.5 GeV, used in §. 

Gross and Surya made a static approximation of all the t- and u- channel graphs. Their 
approximation leads to separability of the potential and hence the complexity of solving 
the integral equation numerically is avoided. It also simplifies considerably the subsequent 
photoproduction analysis, since the meson- and isobar-exchange currents have a form of 
contact terms. On the other hand, the spin-3/2 (and higher) waves, such as P33 and D13, 
receive contributions only from the s-channel exchanges of corresponding resonances, which 
in particular leads to overestimates of the resonance coupling parameters. Also, the static 
(zero-range) approximation of the u-exchange potential can be justified only at low energy, 
because in principle the range of such a potential rapidly increases from 1 / v^2mm^ (where 
m is the exchanged particle mass) at threshold till 1 / m,r at high energy. We observed that 
already at 100 MeV pion kinetic energy the solution of the integral equation for the static 
or exact w-exchange potential differ significantly. In contrast, the range of the t-exchange 
potential is determined only by the exchange mass, and if that mass is heavy enough the 
static approximation may be applicable. Another difference comes from the fact that Gross 
and Surya use the admixture of pseudoscalar and pseudovector coupling for the vrA^A^ vertex. 
Consequent differences, motivated by the consistency with the soft-pion limit, appear in the 
form of the p and a couplings. 

Lahiff and Afnan do not include resonances beyond the A(1232) but apart from that 
they use an interaction, which is very similar to ours. They have as well compared the con- 
ventional versus gauge-invariant vriVA coupling. However, in contrast to us, they preferred 
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the former one, particularly because the spin-1/2 background adds a significant attraction in 
the Pii channel, which helps to simulate partially the Roper-resonance behavior. We include 
this resonance explicitly and hence the spin-1/2 background is not necessary for fitting Pn 
phase shift. 

Including the A and Di^ resonance in the Lagrangian via the relativistic spin-3/2 fields, 
we have used couplings which do not involve the spin-1/2 sector of the Rarita-Schwinger 
theory, and therefore no "spin-1/2 backgrounds" associated with these particles appear in 
the S'-matrix. We do not need these backgrounds to obtain a proper description of the 
data. On a simple example (see Tables ^ and [rV| ) we have seen that, even when these 
backgrounds may sometimes seem relevant for the description, their role can be taken by 
other mechanisms, hopefully with more sensible physical interpretation. 

Although the presented model improves in some aspects previous relativistic analyses 
based on the potential approach, the difficulty of controlling the chiral symmetry constraint 
is present here as well. We do begin with a Lagrangian, and thus the driving force, consistent 
with chiral symmetry, but this consistency can in principle be spoiled by rescatterings, 
particularly because of the loss of crossing symmetry. However, numerical checks (see, e.g., 
Ref. 0) indicate that the amount of violation of the soft-pion theorems is usually negligible 
in the potential modeling. Especially if to take into account that the prime objection of such 
models is to describe the nN physics at intermediate energies where unitarity aspects take 
the leading role. It would, of course, be anyhow important to build the chiral constraint 
more precisely into the nonperturbative vrA^ models. 

In conclusion, we have obtained a description of the nN force in a relativistic dynamical 
model based on the covariant equal-time (quasipotential) equation for a hadron-exchange 
potential. The good quality of our fits in the region up to 600 MeV pion kinetic energy 
suggests the used force and the relativistic approach may be considered reasonable, even 
though this model still lacks inelastic mechanisms which can become important in Su and 
Pii channels above 400 MeV. It is clearly of interest to study the model predictions in other 
processes, such as pion photoproduction and Compton scattering in the vrA^ system. 
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APPENDIX A: CONVENTIONS 

Metric g^^ = diag(+l, —1, —1, —1); Levy-Cevita symbols: = 1, e'^''^^ = e''^^. 
Pauli spinors: in the direction 9, ip are given by 

Xa(^,^)= E 4'We^^'-''^"XA'(0), 

A'=-l/2 

where Xi/2(0) — (l^O); X-i/2(0) — (0; 1); ^^"^ '^aa'(^) ^'^^ Wigner li-functions. 
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Positive- and negative-energy helicity spinors: 



where N± = y {Ep ± m)/2Ep, Ep = ^/mP~+~p^, and p, 6, ip are the spherical coordinates of 
p. The hehcity spinors satisfy the following orthogonality and completeness conditions: 

p=± A 



APPENDIX B: PARTIAL- WAVE OFF-SHELL ttN AMPLITUDES 

The partial- wave reduction is done in the CMS, where the total four-momentum P = 
{PqjO). Using the orthonormality and completeness of the nucleon helicity spinors (see 
Appendix A), we write the BS equation for the helicity amplitudes: 



p"\" 



where we have assumed that the vrA^ propagator is diagonal in the helicity basis, i.e., 
which is true in the CMS. Eq. ( |B1[ ) yields (omitting the momenta) 



oo oo 



rpJp'p _ jrJp'p I ^ 



-OO P ^ 



(Bl) 



(B2) 



(B3) 



where 



x(,1{p\p-Po) = E {J + ^)Xifl\p',,v\po.V.P^)diM. 

J=l/2 

^/'<'(p'o,p',po,p;Po) = C d{cose)x(,{{p\p-p,)di,^{e), 



(B4) 
(B5) 



with X = y or T. We have chosen the 3-vectors p and to lie in the XZ-plane (hence 
Lf = ip' = Q), and 9 is the center-of-mass scattering angle. 

Parity conservation infers the following symmetry for the partial-wave helicity ampli- 
tudes: 



T^xip'o^ V\PO: P; ^o) = P pT%'',-x{Po: V\Po, V] Pq)- 
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(B6) 



These relations reduce the number of independent amphtudes from sixteen to eight. Time- 
reversal invariance imphes 



Tw\{p'o^ v\ Po, P; ^o) = 4A'A T^x'^iPo, p', Po, P; ^o) 



Jp'pf I I 



(B7) 



which obviously does not give any new relations. 

It is convenient to introduce the partial-wave state with definite parity r in terms of the 



partial-wave helicity state |51,52 



|J,p, A) +rp |J,p, -A) 
V2 



[B8) 



The BS equation for the parity-conserving amplitudes takes the form 

oo oo 



(B9) 



-oo 



where r denotes parity. Note the simplification: Eq. ( p9|) has two less coupled channels than 
Eq. (P5D. The partial- wave decomposition of the ET equation (|TU]) proceeds in the same 
way. In addition we use in this case the fact that in the CMS V and T are independent of 
go- The resulting equations are given by Eq. (|T^). 

The amplitudes (0) are related in a simple way to the parity-conserving partial-wave 
amplitudes of Jacob and Wick Starting from representation Eq. (|T3p, we first reduce 
the spinor algebra to the subspace of Pauli spinors Xa(^) and use 



With the aid of Eqs. ( p5|) , (p^ ) and the identity 

d\(li0)di,,{9) = i (p_^_i(cos^) +4AA'P^^i(cose)) 



(BIO) 



(Bll) 



we find the sought expression for the partial-wave amplitude T^^'p with definite parity r in 
terms of T±: 



T^p'p = 1 d(cos 9) \TP'p{e) P, 1 (cos 9) Rp'p + T^^iO) P,_^ i (cos 9) R^p 

I J-l L -'~2 -^+2 



(B12) 



where the dependence on external momenta is omitted for brevity and only the dependence 
on the scattering angle is exhibited. The factors 



pp'NpN' r = + 



2Er, 



N' 



Epi ± rrtN 
lEpi 



(B13) 
(B14) 



arise from the helicity spinors. 
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Now we can write down the relations due to the charge conjugation symmetry. It relates 
the amplitudes of positive energy with those of negative energy and opposite parity: 



T/('''^)(p'o,p',Po,p;Po)=T, 



J(-p'-p), 



(B15) 



Needless to say, these relations will only hold in the quasipotential formulations which satisfy 
the charge conjugation symmetry. 

The isospin decomposition of the vrA^ amplitude is carried out as follows. If we denote 
Xn and 0" as the isospin states of the nucleon and pion, respectively, then 



I /a I 



= r.x^N [i^aT, + (5^, _ i^^r,) ^^^l^ (B16) 
where Ta (a = 1, 2, 3) are the isospin Pauli matrices, satisfying ^[r^, T}\ = iSabcTc- Evidently, 



(BIT) 



APPENDIX C: THE OFF-SHELL ttN POTENTIAL 

According to Eq. (|ID the off-shell vrA^ amplitude can completely be specified by the scalar 
2x2 matrices A and B. Here we give the expressions for these matrices corresponding to 
the tree-diagram potential in Fig. H. Also the isospin factors, J-'{I), are given. 

1. Baryon-exchange graphs 



For the isospin- 1/2 baryon: J-'{^) = —1, ^(|) = 2. 
For the isospin-3/2 baryon: J^{^) = 2, ^(|) = |. 

(a) The u-channel exchange of a baryon with spin ^, mass m and parity rj, using the (pseudo- 
)scalar vs (pseudo-)vector admixture coupling, cf. Ref. specified by parameter A [A = 
corresponds to pure (pseudo-)vector and is used in the text]: 



A' 



y-wNB 



Att U 



A 



m r]\ A(l — A) 



+ 



Am? 



B' 



y-nNB 



7] 

m{u — p'^ 



-2riu + tjIp"^ + p'"^) — m 
2m \^ -m 2r] 

/2n r]{p''+p'^) 



-m 



All u 



A 







+ 



A(l - A) / 2m -7] 
2m \-r]0 



+ 



Am? 



( —rju 


m \ 


\ m 





HI)- 



(b) The A exchange using the conventional coupling, Eq. (0): 



(CI) 



(C2) 
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j^I _ fnNA 



47rm^ u — m? 



+ 







2>m? 
2azm 



+ m{\ — 2az) 



p' -Pu- P22 - p'^ \ 

\p-Pu- P22-P^ ) 

P22 - {p'^p + pV) • Pu/u A 

I - {p + p') -pu/u j 

' {p'\ + p^A;') . - (A; + k') ■ Pu - 2P22 

{k + k') ■ Pu/u 



^ f k'-pu + 2P22 

M A; • p„ + 2P22 



1 



- (1 + ayP22 ( 1 J ~ + 2a^)P, 
7 /, 
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1 




P, 



+ 



47rm^ M — I ^2 
2 / -P22 \ 



3m^ 
+ m( 



1 



p'-pu/u 
P ■ Pu/u 



2azm 



{k'+k) ■PU+2P22 




+ (l + afP22 



1 




where 



a = —z — 



2' 



P22 = {k' ■ Pu) (k ■ Pu)/u, P22 = (p ■ Pu) {p ■ Pu)/u, 
Pu = p - k' = p' - k, u = pI, 

and A 3, B 3 are the contributions of the spin-3/2 projection operator. 

^2 ^2 

For half-integer spin the contributions of the spin-j projection operator read 
1 



(2j)! 



1 or^+i^""^ 



m[P22 - {p''^P + pV) ■ Pu/ u] p' -Pu- P22 - P' 



+ 



/2 



P-Pu- P22- P 



Bpj = 



+ 



(2j)!! 

P22 mp'W.U, • 
mp-pu/u 1 / j-2' 







] i+i^^"'' 



where P/(x) is the first derivative of the Legendre polynomial, and 



25 



r _ , jk-Pu) 7, _ w {k' ■ pu) _ k-k' 

K — K K — K Pu-i Xu - 



U 



U 



(c) The A exchange using the gauge- invariant couphng (|25|), according to Eq. 



A' 



tNA 



u 



AirmT^m u — rri 
fnNA "^^ 



— rr)2 P 



2 P 



(C5) 



nves 



(C6) 
(C7) 



2. Meson-exchange graphs 



For the isoscalar meson: J^{^) = ^(|) = 2. 
For the isovector meson: J^(^) = 2, = — 1. 

(a) p-meson exchange: 



A 



I gpTTirSpNN 



^{e + k'^-2s) -1 



47r rn?p — t 

9pmT9pNN 1 

Att — t 



2mjv 



-1 



(b) (T-meson exchange: 



/ QcjttttQctNN 



A 



B^ = 0. 



-Xk^ + k'^) 



mi — t 



1 




HI), 



(C8) 
(C9) 



(CIO) 

(Cll) 



APPENDIX D: ttN AMPLITUDES FOR HIGHER-SPIN BARYON EXCHANGE 



The use of GI couphngs of higher-spin baryons |31|] allows us to treat exchanges of 
baryons with any spin in a straightforward way. Taking the point of view that consistent 
ttNN* couplings (where by N* we denote the spin-j baryon, j > 3/2) are those invariant 
under the appropriate gauge transformations of the N* field, we can write down the following 
ansatz for the nN scattering amplitude through a spin-j baryon tree-level exchange: 



M{k',k;P) 



n 



' (2/ + 1)!! ^-m 
P^^^^-k.k')+}im-k-k'] 



P^\Jk^k/^ 



(Dl) 



where / = j ~\i Pi{x) is the first derivative of the Legendre polynomial, P is the momentum 
of the exchanged baryon, and 
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p2 



kfj, — k^/yk 



(D2) 



This amphtude is actually just the spin-j projection operator contracted with the external 
pion momenta and multiplied by g'^p'^\'][) — m)~^ . Because of the projection operator, the 
amplitude contains only the spin-j contributions, as we will now demonstrate. 

In the CMS P = (W, 0), hence k = (0, k) = (0, —p). The s-channel helicity amplitude is 
then written as, 



M'y'x = 4il^)M{k',k;p)ui{p) = C,u(,ip' 



(D3) 



where Xs = cosd, 9 is the cm. scattering angle, s = W'^, and the constant factors are 
absorbed in Cj. Using the Dirac equation, {^p)u1{p) = (p-Ep 7o — ttT'n)u1{p), we obtain 



-H^7o + m p'Epi^Q - TUN pEp'jo - tun 



x'x — ""a' (p' 
+ 



s — 



P 



ipp'sYuiip). 



(D4) 



The contribution of this amplitude to the partial wave with total spin J is given by the 
corresponding partial-wave amplitude. The latter can be obtained by the procedure of 
Appendix B, for that we only need to know the following angular integrals: 



dxPl(x)PL(x) 



1, L = l-2n-l 
0, otherwise 



(D5) 



where n should be an integer between and (Z — l)/2. The resulting parity-conserving 
partial-wave amplitude is then found to be 



J = j, 



MP'P^ = Cj (pp's)' X \ {R^y, + M'_,M.r R'r'). J = J - 1 - 2n, 

{Rpy + ACA/; R^yr) , J = J - 2 - 2n, 



(D6) 



where A/± = {pEp ± m^)/p, S^'^'^ = {±W — m)^^, and factors R are defined in Eq. ( |B13D . 
Using the explicit form of R, we can see that the lower partial-wave contributions vanish 
exactly. Thus, amplitude Eq. (pi|) has only the highest-spin contribution: 



Mfp^ = 6,jCj{pp'sys^''^ Rpy 



(D7) 
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FIGURES 
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FIG. 1. Diagrammatic form of a relativistic vrA^ scattering equation. 
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FIG. 2. The tree-level vrA^ potential. 
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FIG. 3. Cutoff form factors as a function of the off-shell 3-momentum squared. 
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FIG. 4. Cutoff form factors as the function of the pion kinetic energy in the lab. 
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FIG. 5. Rewriting the equation for the case when the potential can be presented as the sum 
of separable and non-separable potentials. 
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(q/mj 

FIG. 6. Effective range plots for the nucleon and p-meson exchange contribution. Dashed hnes: 
Kp = (Weinberg- Tomozawa contact term). Solid lines: Kp = 3.7 (vector-meson dominance). The 
data points are extracted from the KH80 (stars) and SM95 (dots) partial-wave analyses. 
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FIG. 7. Effective range plots for the sum of the N, A and p exchanges with g^^^/Ai^ = 13.6, 
g'p/47r = 3, Kp = 3.7, /^ata/^'^ = 0.36. Dashed hnes: Peccei choice; dashed-dotted lines: NEK 
choice; solid lines: gauge-invariant vrA^A coupling. Results for ^n, S'si, and P33 waves are shown 
in bold lines. Data points are same as in Fig. |6[ 
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FIG. 8. The description of the S- and P- and D-wave ttN phase-shifts up to 600 MeV. 
Data points are from the SM95 (triangles) and KH80 (stars) partial-wave analysis. Solid lines 
represent the model solution. Dotted lines represent the calculation where the principal value part 
of the rescattering integrals is switched off (i.e., the K-matrix approximation with the same set of 
parameters). Dashed line for the Sn shows the calculation when the pole contribution of the Sn 
resonance is switched off. 
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FIG. 9. The description of vrA^ phase-shifts up to 350 MeV pion lab kinetic energy. Sohd hnes 
represent the fuU model calculation. Legend for the data points is the same as in Fig. pi 
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FIG. 10. Description of the effective range plots for the S- and P-waves. Solid hnes represent 
model calculation. Dotted hnes are the ii'-matrix approximation with the same set parameters. 
Dashed-dotted line represent the full calculation with Kp = 3.7. 
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FIG. 12. The real and imaginary parts of the self-energies 
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FIG. 13. The energy dependence of the A'^A'^* mixing angle. The numbers near the curve stand 
for the corresponding value of the pion lab kinetic energy in GeV. 
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FIG. 14. The dynamical ttNN, ttNA and nNDi^ form factors, as a function of the off-shell 
3-momentum squared of the ttN state. The energy of the baryon is fixed at the ttN threshold. 
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FIG. 15. The dynamical half-off-shell form factors for the case when the pion and nucleon are 
on-shell, as a function of the pion lab kinetic energy. 
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TABLES 



hi2j 


Np (WT) Np (VMD) 


Peccei 


NEK NpA (WT) NpA (VMD) 


KH80 m 


KA86 [0 


SM95 1 


32| 


Su 


0.171 


0.171 


0.144 


0.171 


0.171 


0.171 


0.173 


0.175 


0.172 


S31 


-0.100 


-0.100 


-0.150 


-0.100 


-0.100 


-0.100 


-0.101 


-0.100 


-0.097 


Pll 


-0.185 


-0.145 


-0.077 


-0.076 


-0.127 


-0.092 


-0.081 


-0.078 


-0.068 


P31 


-0.059 


-0.079 


-0.064 


-0.085 


-0.043 


-0.066 


-0.045 


-0.043 


-0.040 


Pl3 


-0.047 


-0.064 


-0.049 


-0.069 


-0.032 


-0.050 


-0.030 


-0.030 


-0.021 


P33 


0.101 


0.109 


0.194 


0.197 


0.224 


0.222 


0.214 


0.215 


0.209 




177 


129 


35 


31 


26 


15 


2.1 


1.6 




X 



TABLE L The S- and P-wave vrA^ scattering lengths in different models are compared to the 
partial-wave analyses. The represent the chisquare value per point with respect to the SM95 
analysis (the analysis error is taken to be the same for all points and is equal to 0.005). Parameters 
corresponding to different calculations are given in Table 



parameters 


Np (WT) Np (VMD) 


Peccei 


NEK 


iVpA (WT) 


iVpA (VMD) 




13.8 13.8 


13.8 


13.8 


13.8 


13.8 




5piVAf5p7r7r/47r 


3.15 3.15 


3.4 


3.15 


3.15 


3.15 






0.0 3.7 


3.7 


0.0 


0.0 


3.7 






X X 


0.31 


0.36 


0.6 


0.55 






X X 


-0.25 


0.5 


X 


X 




TABLE IL 


Parameters corresponding to calculations of Table |. 


Values in bold 


were varied to 


give a best fit of the scattering lengths. 












hi2j 


NpAa (A) 




iVpAcr (B) 


SM95 


321 


Sii 


0.170 




0.170 




0.172 


S31 


-0.100 




-0.100 




-0.097 


Pll 


-0.068 




-0.069 




-0.068 


P31 


-0.039 




-0.039 




-0.040 


Pl3 


-0.024 




-0.023 




-0.021 


P33 


0.209 




0.209 




0.209 




0.16 




0.13 


X 



TABLE III. Scattering lengths resulting from the sum of N ^ A, p, and a tree-level exchanges. 



Parameters corresponding to different models are given in Table IV. 
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parameters 



iVpAfj (A) 



iV/>AfT (B) 



91nn/^'^ 

gpNNgpTTTr/^'^ 
gaNNgaTTTrmN/ (STTm^) 

Crr 



13.8 
3.15 
0.0 
0.28 
0.5 
2.1 
1.0 



13.8 
3.15 
4.6 
0.30 

X 

2.0 

1.0 



TABLE IV. Parameters corresponding to the calculation in Table |II|. Values in bold were 
established by the fit. Model A uses the conventional vriVA coupling, while model B uses the 
"gauge-invariant" coupling. 



field 



coupling constants 



An 



masses [GeV] 



N 

TT 

N* 
A 

P 
a 

Sii 



13.8 (/2^^=0.0757) 



9ttNN* _ 2 7 

0.43 



47r 

J-irNA 



nND. _ 0.6 



47r 
QpNNgpTT 
An 

aoNNSlaj^nif^ = 1.0, = 0.75 
2m,r 



2.85, Kp = 1.8 



An 
An 



0.6 



rriN = 0.9383, yl^v = 1-8 
= 0.139, = 1.16 

niN* = 1.438, Am* = Am 
rriA = 1.252, Aa = 1.5 
niD = 1.525, Ad = 1.7 
TUp = 0.77, Ap = 1.9 
m^j = 0.55, Afj = 1.1 

msi, = 1.555, /l5n = 2.0 



TABLE V. The model parameters. The values given in bold were varied during the fit. 



field 



coupling constants 



An 



masses [GeV] 



N 

TT 

N* 
A 

^13 
P 

a 



13.8 (4V=0.0757) 



= 2.6 
0.36 



2 

S-kNN* 
An 

An 

JtvND 
An 

gpNN9pTTw 
An 



0.8 



3.1, Kp = 2.7 



An 



2mTT 



1.3, = 1.0 



mjv = 0.9383, Am = 2.2 
= 0.139 

niM* = 1.47, Am* = Am 

ruA = 1.232, tIa = 2.6 

ruD = 1.525, Ad = 1.7 
rup = 0.77, Ap = 1.3 
= 0.55, Afj = 1.9 



TABLE VI. The parameters of the ii'-matrix calculation, see Fig. 11. 
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field 


Zi 


Z-i 


mn [GeVl 


N 


0.63 


0.72 


1.090 


N* 


0.63 


0.95 


0.960 


A 


0.68 


1.01 


1.390 




0.93 


1.03 


1.512 


Su 


0.97 


1.09 


1.555 



TABLE VII. The renormalization constants obtained in the model. 
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